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MgB 2 under pressure: phonon calculations, Raman spectroscopy, 

and optical reflectance 
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Max-Planck-Institut fur Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 

(received 21 May 2001) 

The effect of pressure on optical phonon frequencies of MgB2 has been calculated using the frozen- 
phonon approach based on a pseudopotential method. Griineisen parameters of the harmonic mode 
frequencies are reported for the high-frequency zone-center E2 g and Bi g and the zone-boundary 
E2u and B2u modes at A. Anharmonic effects of phonon frequencies and the implications of the 
calculated phonon frequency shifts for the pressure dependence of the superconducting transition 
temperature of MgB2 are discussed. Also reported are Raman and optical reflectance spectra of 
MgB2 measured at high pressures. The experimental observations in combination with calculated 
results indicate that broad spectral features we observed in the Raman spectra at frequencies between 
500 and 900 cm -1 cannot be attributed to first-order scattering by zone-center modes, but originate 
in part from a chemical species other than MgB2 at the sample surface and in part from a maximum 
in the MgB2 phonon density of states. Low-temperature Raman spectra taken at ambient pressure 
showed increased scattering intensity in the region below 300 cm -1 . 



PACS: 63.20-e, 74.25.Kc, 78.20.Ci, 78.30 .-j, 62.50+p 



I. INTRODUCTION 

Magnesium diboride, MgB2, was recently discovered to 
exhibit a superconducting transition at 39 K Q , which is 
by far the highest T c for a binary compound. The light 
atomic masses in MgB2 enhance the phonon frequencies 
which set the scale for T c in BCS theory. A variety of ex- 
perimental observations (e.g., isotope effects and scan- 
ning tunneling measurements of the superconducting gap 
|^-^|), indeed indicate phonon- mediated (BCS) super- 
conductivity, and theory finds medium or strong electron- 
phonon coupling ^-^|. The coupling is predicted to be 
particularly strong for the zone-center optical phonon of 
E<i g symmetry and the related phonon branch connecting 
to the zone-boundary E<i n mode at the A point (direction 
along the c axis of the AlB2-type structure) . 

The effect of pressure on the superconducting proper- 
ties of MgB 2 has been studied by several groups. |3~0|— p~^] 
While there is some scatter in the experimental pressure 
coefficients all studies show a decrease of T c , at rates 
between 0.7 and 2.0 K/GPa. Lattice parameters un- 
der pressure have been determined by neutron and x- 
ray diffraction |]l6|-[l9|. An ab initio calculation of the 
pressure- volume relation and the optimized c/a axial ra- 
tio under pressure ]2C[| is in good agreement with the 
experimental data. It was concluded that the decrease in 
T c is not driven by a change in the electronic density of 
states near Ep but by an increase in phonon frequencies, 
specifically of the strongly coupled E^ g mode . 

The frequency of the E2 g (T) phonon is expected be- 
tween 470 and 665 cm -1 , according to recent calcula- 
tions P,p2Pj2^,p| . So far, this phonon mode has not 
been identified unambiguously by neutron inelastic scat- 
tering HHHI or in Raman experiments ffl|3|^^Q. 
A Raman spectrum by Bohnen et al. [23] indicates a 
broad asymmetric feature near 580 cm , while other 



authors (lJJ^,^! report an even broader feature (260- 
300 cm -1 width) centered near 620 cm -1 . Contradicto- 
rily, compliance with the selection rules for a Ei g mode 
has been reported |27| as well as complete breakdown of 
the same, attributed to resonance effects p(| . The 620- 
cm _1 feature was observed to increase in frequency under 
pressure, with an unusually large Griineisen parameter 
Jl9| . It was generally assigned to first-order zone-center 
phonon scattering, highly broadened by strong electron- 
phonon coupling or anharmonic effects, but other inter- 
pretations were also considered. 

In this paper we present a theoretical investigation of 
the pressure dependence of phonon frequencies of MgB2. 
As no phonon has been, so far, unambiguously identified 
experimentally, the theoretical methods gain in impor- 
tance - in particular the ab initio approaches, such as the 
frozen phonon method which is employed in this work. 
Compared to linear response theory, the frozen phonon 
method has the disadvantage that it is not very efficient 
in generating an overall picture of the phonon dispersions 
throughout the Brillouin zone, but is restricted to calcu- 
lation of a few high-symmetry modes. One of the advan- 
tages of the frozen phonon method, however, is that it 
yields a clear picture of anharmonic effects, which in ear- 
lier calculations || were found to be particularly large 
for the atomic displacements corresponding to the Ei g 
mode of MgB2. The immediate interest in the shifts of 
phonon frequencies is in the context of interpreting the 
pressure dependence of T c . Furthermore, the calculated 
pressure dependences of phonon frequencies can in prin- 
ciple be useful for the identification of spectral features 
in Raman spectra. 

We also report here Raman and optical reflectance 
measurements of MgB2 as a function of pressure at room 
temperature. Based on the pressure dependence of Ra- 
man spectra and the optical reflectance and by taking 
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into account the results of the phonon frequency calcu- 
lations we argue that the observed Raman features are 
not due to zone-center first-order Raman scattering from 
MgB2 . The interpretation we suggest for our Raman re- 
sults considers a superposition of a density-of-states peak 
of MgB2 and scattering by a different chemical species 
present at the surface of as-grown samples. We also take 
a look at zero-pressure Raman spectra measured as a 
function of temperature in order to check whether it is 
possible to observe a pair-breaking excitation below T c . 

The paper is organized as follows: the symmetry anal- 
ysis of all phonon modes at the T and A point is presented 
in Section EJ. Details of the theoretical method and the 
calculated results are presented and discussed in Section 
Experimental details and results are given in Section 
followed by conclusions in Section M. 
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II. SYMMETRY ANALYSIS 

The frozen phonon calculations described below re- 
quire knowledge of the phonon displacement patterns. 
We thus start with the symmetry analysis of the vibra- 
tional spectra. The point group of the AIB2 structure is 
D^fr. A factor-group analysis yields the decomposition of 
the coordinate representation 



D 



is 



E 2g + 2A 2u + 2E lu at T 



and 



Aig + E\ g + A 2u + B 2u + E\u + E 2u at A 



(1) 



(2) 



with, at the point T, one A 2u and one E\ u referring to the 
rigid translations (w = 0). Constructing then projection 
operators we find the bases spanning the irreducible rep- 
resentations in Eqs. 1 and 2 - i.e., the respective phonon 
displacement patterns at T and A, respectively. They 
are shown in Figs, [l] and |^, with the modes arranged 
(columnwise) in the order of decreasing MgB2 eigenfre- 
quencies (which are obtained later, in Section II. C). At 
high pressures this ordering may, of course, be altered. 
The patterns in Fig. || are arranged so as to be the "end- 
points" of those T-A branches of the phonon dispersion 
that start with the corresponding pattern in Fig. [l| This 
means that e.g. the three patterns in the last column of 
Fig. [2] are the end-points of the acoustic branches. 

The patterns at A alternate the +u and — u displace- 
ments from one layer to another. In addition, for the two 
highest frequency modes we notice that (to within the ± 
alternation) the patterns of a layer are the same as at T. 



III. CALCULATIONS 



A. Theoretical Method 

The frozen phonon approach starts from the ab ini- 
tio evaluation of the total energy E tot of the solid with 
frozen-in atomic displacements. The energy is evalu- 
ated using the density functional theory (DFT) within 
the generalized gradient approximation (GGA) p3j. We 
work in a plane- wave basis and use pseudopotentials. 
For the actual calculations we employed the VASP codes 
p9| p2[ , and the ultra-soft Vanderbilt-type pseudopo- 
tentials |33| were supplied by Kresse and Hafner Q. 
The pseudopotential we chose for Mg treats the semi- 
core states explicitely, i.e., Mg is represented with 8 va- 
lence electrons 2p 6 3s 2 , and the non-linear core correction 
(NLCC) [^H is applied as well. The calculations are car- 
ried out with the plane- wave cutoff energy of 33.6 Ry, 
and the Brillouin zone sampling is based on T-centered 
uniform meshes 16 x 16 x 14 or 16 x 16 x 8 - the latter 
one applied in the calculations with the doubled unit cell 
for frozen phonons at the Brillouin zone boundary. As 
the system is metallic, the fc-space integration with the 
incompletely filled orbitals uses the tetrahedron method 
f{6"f with BlochPs corrections |37) . The above meshes di- 
vide the Brillouin zone into 21504 or 12288 tetrahcdra; 
the number of irreducible tetrahedra is lower by a fac- 
tor 2 to 6, depending on the symmetry of the frozen-in 
phonon. 



B. Structural Properties 

For eleven volumes between 25.0 A 3 and 29.6 A 3 , we 
calculated E tot (V) with different c/a values while keeping 
the volume V constant and, for every chosen V, found the 
optimized structure, i.e. the corresponding a(V), c(V), 
c/a(V), and E(V). The latter data were fitted by the 
Murnaghan relation for E(V) |3^] which provided the 
static equilibrium J39| as well as the bulk modulus Bq and 
its pressure derivative B'. The results are summarized 
in Table 0. The calculated equilibrium volume (Vb = 
28.663 A 3 ) compares well with the experimental value 
V = 28.917(1) A 3 (Ref. pfl, at 37 K) meaning that we 
are -0.9% off in Vq. The individual lattice parameters are 
-0.6% off the experiment for ao, +0.2% for c , and +0.8 
% for co/do- The variation of c/a under pressure has 
been translated into a quadratic function of (1 — V/Vq) 
and its coefficients are given in Table | as well. The 
present results are also in good agreement with those of 
a previous calculation |^] (see Table |). 

A limited testing of different pseudopotentials revealed 
that treating the magnesium 2p 6 electrons as valence 
states is not essential for getting the correct equilibrium, 
but using the GGA (rather than the LDA) is [|cj . With 
the value of c/a = 1.14 the MgB2 is well inside the 
stability region of the AIB2 structures H|. A rather 
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high pressure would probably be required to induce a 
phase transition, possibly towards the UHg2-type struc- 
ture (isostructural to A1B 2 , but c/a less than ~ 0.8) or a 
variant thereof, as is observed in other Allotype inter- 
metallic compounds p2[ . 

The subsequent calculations of phonon frequencies are 
performed at the calculated equilibrium, and at volumes 
related to pressure by the Murnaghan relation 



P{V) = § 



(v/v )- 



(3) 



which uses the calculated Vq, Bq, B' quoted in Table |]. 



C. Phonon frequencies 

In this work we are mainly concerned with the two 
highest frequency modes Bi g (T) and E2 g (T) and their 
A-point counterparts B 2u (A) and E 2u (A). For compari- 
son with previous works we calculated the A 2u (T) and 
E lu (T) as well. The branch E 2g (T) - E 2u {A) in the 
phonon dispersion is the one that exhibits the strongest 
electron-phonon coupling ||, and the Bi g (T) is the z- 
displacement analogue of the E2 g (T). 

We determine the vibrational eigenfrequencies using 
the frozen phonon method [ ff3|] , while paying particular 
attention to obtaining the harmonic part of the ab initio 
calculated energy differences that, inevitably, have to be 
probed by finite displacements u. Thus for each pattern, 
atoms are given 5 different displacements ranging from 
u/a = 0.010 to 0.035, and the (six) calculated energy 
values E{u) are fitted with a quartic polynomial yielding 



AE tot {u) = a 2 (u/a) 2 + a 3 (u/a) 3 + a 4 (w/a) 4 (4) 

As an example we show in Fig. ^ the variation of E tot (u) 
for the E2 9 (r) mode 44 1. The harmonic part of AE to t is 
then retained for determination of the phonon energy: 



uj 1 

— M K u 2 K = E harm {u ^ 0) - E harm (u = 0) (5) 

For the four phonons listed above the procedure is re- 
peated at various volumes. Although straightforward, 
this procedure calls for a few explanations. 

The frozen-in phonon displacements of Fig. ^ conserve 
the translational symmetry of the elementary unit cell 
defining the MgB2 structure. With the zone boundary 
modes given in Fig. || the total energy calculations as 
a function of displacement u are to be done on a unit 
cell obtained by doubling the elementary cell (viz. along 
the [001] direction), which then comprises two formula 
units of MgB2. The frozen-in displacements also lower 
the rotational symmetry of the system: starting from 



the D 6 h, which is the point- group of the MgB 2 lattice, we 
end with, respectively, D 3 d, D 2 h, C 2 h when the Bi g (T), 
^2g(r)-(a,b) displacements are imposed, and D 3 h, D 2 h, 
D 2 h with the B 2u (A), E 2u (A)-(&,b) displacements, re- 
spectively. 

As the displacements u are finite, a larger or smaller 
degree of anharmonicity will always be present, but some 
restrictions on the latter can be read out from symmetry. 
With all patterns at the zone boundary, replacing +u by 
— u has the same effect as shifting the phonon distorted 
structure by the lattice constant c, and it leaves the E{u) 
invariant. Consequently, the cubic (and, generally, odd- 
power-) anharmonicity is forbidden in all zone-boundary 
modes. Inspecting, in turn, the T-displacement patterns 
we realize that for nearly all modes the +u <-> — u sub- 
stitution is equivalent to applying a mirror reflection or 
C 2 rotation - an operation which, in spite of the low- 
ered symmetry, is still present in the point groups of the 
phonon-modulated structures - and which thus leaves in- 
variant the structure and its energy E(u). This forbids 
the cubic term in the expansion of E(u) as well. Among 
the 9 modes in Fig. [l] the E 2g (T)-(a) pattern is the only 
one in which all anharmonic terms are allowed and in 
which the cubic anharmonicity is to be expected. 

For evaluation of phonon frequencies from AE to t(u) we 
used the atomic masses of B and Mg that correspond to 
the natural isotope distribution: 10.8 and 24.3 a.m.u.; 
this is the usual approach in the spirit of the virtual 
crystal approximation. We may keep in mind that for 
a 10 B-enriched sample the frequencies of all the modes 
involving boron would be higher by 4%. 

Table summarizes the results of our frozen phonon 
calculations at different volumes: the variations with vol- 
ume or pressure are described by the mode Griineisen pa- 
rameter 70 = — [d In uj/d In V] Vq or the coefficients a and 
6 of a quadratic fit (see the caption) . Figure [| shows the 
calculated phonon frequencies as a function of volume. 

We estimate that the calculated harmonic frequencies 
are uncertain to ~l-3% (depending on the mode); this is 
the "instrumental precision" of the determination of to- 
tal energies and the uncertainty following from the error 
analysis in fitting the variations E(u) by polynomials. 



D. Anharmonicity and the E2 g (T) mode 

We note that the five displacements u/a quoted above 
range from one half to approximately the double of the 
r.m.s. displacement [m], and one expects vibrations in 
both harmonic and anharmonic regimes. The effort of re- 
moving from AEtot the anharmonic contributions yields, 
in turn, an insight into the anharmonicities. 

The phonons exhibiting the strongest electron-phonon 
coupling are those of the E 2g (T) — E 2u (A) branch H. 
The E 2g (r) has already received considerable attention 
2^J3J23||| . Being doubly degenerate, it is materialized 
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by two distinct displacement patterns which are labeled 
(a) and (b) in Fig. [|. They both should lead to the same 
eigenfrequency (the same coefficient a 2 in Eq. ||) and, in 
this sense, they are equivalent - as long as we stay within 
the harmonic approximation, i.e., within the standard 
group-theory reasonings. Once we step out of the har- 
monic approximation, the displacements (a) and (b) will 
become two unrelated patterns which are described by 
different anharmonic coefficients. This gets confirmed by 
numerical calculations and we obtained a 2 = 107 eV/cell 
with both (a) and (b) patterns and, somewhat unexpect- 
edly, also 04 = 45800 eV/cell equal in both patterns. The 
two modes are nevertheless distinct in the cubic term, 
which is absent in the (b) pattern (03 = 0) and given by 
a 3 = -f644 eV/cell in the (a)-pattern. 

It is not difficult to understand why the mode -E 2g (r)- 
(a) exhibits a strong cubic anharmonicity: giving the 
boron atoms a +u or — u displacement, the main con- 
tribution to AE tot comes from stretching or compressing 
the B-B bonds - which, obviously, has very different en- 
ergetic cost. In the calculations, the cubic anharmonicity 
of the mode is dealt with easily: after obtaining the AE to t 
with the displacements +u and — u we take an average, 
and are left with AE to t carrying quartic anharmonicity 
only. Using the pattern (b) instead of (a) is technically 
somewhat simpler, because the averaging of the contribu- 
tions from the stretched and compressed bonds is "built 
in" within the same displacement pattern. 

A striking feature of the E2 g (T) mode (one which al- 
ready has been noticed by other authors J9|) is its strong 
anharmonicity 04; it is larger by an order of magnitude 
than that of the Bi g (T) mode. This is equally true for 
the corresponding mode £? 2u (A) at the Brillouin zone 
boundary and, apparently, holds for the whole branch 
^2ff(r) — E2u{A). The ratio of the quartic and harmonic 
energies, at any given displacement u, is 04/02 = 4 eV" 1 
for E 2g (T) and 3.5 for E 2u (A). Yildirim et al. Q arrived 
at an even stronger anharmonicity of the E2 g (T) mode: 
di/al « 8 eV" 1 . 

Large anharmonicity is consistent with a strong varia- 
tion of the phonon frequency with volume and the mode 
Griineisen parameters of -E 2g (r) and E2 U (A) turn out to 
be much larger than those of the other modes (see Table 
||) - and, in any case, unusually large: 2.5 and 2.8, re- 
spectively. The corresponding frequency shifts Auj oc jlu 
are approximately the same, which means that the whole 
branch E2 g {T) - E2 U (A) is shifting "rigidly" when com- 
pressed - unlike the branch Bi g (T) — B2 U (A) which be- 
comes more dispersive under pressure. 

It is interesting to note that for the E2 g (T) and i?2« (^4) 
modes the quartic coefficient 04 of Eq. |] (not given in 
Table |n|) turns out to be approximately constant and 
independent of volume, in the entire range of volumes 
shown in Fig. |J this means that, under pressure, the 
quartic energy goes as the inverse fourth power of the 
lattice constant a. 



A well known effect of anharmonicity is in broaden- 
ing the spectral features in phonon spectroscopies and in 
shifting them by some amount A(oj) with respect to the 
ideal 5(lu' — w)-like response of the perfectly harmonic 
system. The intricate details of the variation A(w) are 
beyond the scope of this work. Yildirim et al. || re- 
port that for the E 2g (r)) mode the A(uj) amounts to 
a frequency shift by 17 - 25% above the harmonic fre- 
quency. The estimate is based on their results for a 2 , 04, 
uj(E2 g (T))— 486 cm" 1 and employs two different ways 
of "renormalizing" the harmonic frequency. Taking into 
account that the anharmonic shift, in the lowest order 
of perturbation theory, is proportional to the ratio 04/02 
and to the phonon energy fkv itself [|45f , a simple rescaling 
of the Yildirim et al. results suggests, for our case, a 9 to 
13% upshift of the E2 g (r)-frequency; similarly, the shift 
in E2u{A) can be estimated at 7 to 10% (see the hatched 
areas in Fig. ||). Treating in the same way also the cal- 
culated harmonic phonon frequencies in compressed unit 
cells wc estimate that these shifts reduce the Griineisen 
parameters of the two modes from 2.5 and 2.8 to 70 — 
2.0 - 2.2 and 70 = 2.3 - 2.5, respectively. 



E. Comparison with other calculations 

There have been several calculations of zone-center vi- 
brations [B]|2^|J and phonon dispersion [H|23|,|9]l . The re- 



sults obtained for modes at T are compared in Table III 



The most conspicuous feature of the comparison is the 
large spread of the results for the E2 g (T) mode - the cal- 
culated frequency varies from 470 cm" 1 Q to 665 cm" 1 
p2| - although otherwise good agreement between differ- 
ent authors is found for the other T modes. 

Bohnen et al. [p3| attribute the disagreement in the fre- 
quency of the -E 2g (r) mode to an insufficiently converged 
/c-point sampling in some of the calculations. In our case 
we have ensured that full convergence is reached. 

Large anharmonicity of the -E 2g (r) was proposed as 
a possible cause of the discrepancies by Yildirim et al. 
H . The proper "separation" of the harmonic part in the 
calculated AE to t indeed matters: evaluating the frozen 
phonon energy at u = 0.057 A (=the r.m.s. displacement) 
we obtain, for the i?2 5 (r)-(a) mode, the harmonic term 
in Eq. | of 37.84 meV, the cubic 10.57 meV and the 
quartic 5.33 meV. Billing the quartic contribution into 
the harmonic energy would shift the frequency up by as 
much as 7%, assuming that the cubic term has been cor- 
rectly eliminated (else the frequency would be shifted by 
another ±14%). 

When linear response theory is used for calculation 
of phonon frequencies the results are free of anharmonic 
contributions by construction. Nevertheless the results of 
Refs. [p2|]8U2£| for the E2 g (X) mode show a large spread, 
too (cf. Tab. 
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We note that there is still another factor which can 
be a source of considerable discrepancies. The large 
mode Griineisen parameter implies a high sensitivity of 
the calculated E2 g (T) phonon frequency to the (equilib- 
rium) volume used. Some authors choose to perform 
the phonon calculations at the experimental equilibrium 
Vq XP rather than at the calculated equilibrium V$ h . Both 
choices are legitimate and, as a rule, they do not lead to 
any marked disagreements; the modes with large 7 may 
nevertheless constitute an exception. Our present cal- 
culations refer to the calculated equilibrium Vq , and the 
E2g (r) frequency at the experimental Vq XP would be only 
12 cm -1 lower; the difference however becomes 35 cm -1 
in Ref. [|3| where the authors were careful enough to de- 
termine two sets of frequencies for all T-point phonons, 
viz. one at the theoretical and one at the experimental 
equilibrium. Generally, the existence of two results (one 
at VQ h and one at Vq XP ) is an intrinsic uncertainty un- 
derlying the concept of frozen phonons itself jlt| and, 
in a sense, it can be thought of as a consequence of the 
DFT-LDA. If we estimate that, in a typical DFT-LDA 
calculation, the Vq' 1 is between 1 and 3% under the exper- 
imental value Vq XP , and with the Griineisen parameter of 
the order 7 = 2.5, the two calculations of frequencies of 
E 2g {T) can differ by as much as 2.5 to 7.5% (« 40 cm -1 ). 

Another possible source of discrepancies between 
frozen-phonon results are the errors of statistical nature 
that relate to the details of the fitting of Eq. ^, e.g., the 
choice of the degree of the polynomial and the selection of 
the displacements u for which the -E(w)'s are calculated 
and the polynomial then fitted. Even if the consequent 
differences are relatively small, both choices should be 
clearly stated for the sake of reproducibility. 

The E2 g (T) mode is, apparently, very sensitive to all 
details of the calculation. Besides the convergence with 
the plane-wave cutoff, the fc-point sampling, and the 
technique used for dealing with the incompletely filled 
orbitals (the "metallic sampling" algorithm) - which are 
amenable to quantitative testing - still another source of 
divergences between different calculations is in the choice 
of the pseudopotentials or of the method for determina- 
tion of the total energy. The consequent uncertainties 
are difficult to estimate though. 

Considering the large spread, among different authors, 
in the calculated values of the E2 9 (r) mode frequencies, 
experimentally, one might attempt to discern this mode 
from the other ones by its unusually large value of the 
Griineisen parameter 7 rather than by its frequency. 

F. Mode Griineisen parameters and pressure 
dependence of T c 

The whole optical phonon branch connecting the 
E2 9 (r) and the E2«(A) modes was shown to exhibit 



strong electron-phonon coupling |§,^3|. Within a me- 
dium or strong coupling scenario it is considered to 
largely determine the superconducting properties of 
MgB2. Based on early calculations of the electron- 
phonon coupling constant A « 0.7, an assumed mode 
Griineisen parameter 7 = 1, and a calculated pressure 
dependence of the density of states at the Fermi level of 
d\nN(0)/dP = — 0.31%/GPa, the pressure dependence 
d\nT c /dP of the critical temperature was estimated to 
fall in the range -2.3 to -3.6%/GPa @. The esti- 
mate compared well with experimental data and thus 
supported the notion of BCS superconductivity in MgB 2 . 
The pressure-induced increase in characteristic phonon 
frequency was identified as the dominant cause of the 
decrease of T c under pressure. 

With a refined calculated value of A = 0.87 [||, ujq = 
500 cm -1 , and a typical Coulomb pseudopotential of 
(i* =0.1 we obtain T c> P=0 = 36.7 K based on the McMil- 
lan expression in the form given in Ref. pofl . Following 
the analysis of Ref. |^(| but using an average 7 = 2.3 for 
the E2 g (T)-E2u{A) phonon branch (taking into account 
anharmonicity) and a bulk modulus Bq = 145 GPa as 
calculated above, as well as d\nN/dP = — 0.31%/GPa 
Q, we estimate dT c /dP « -1.8 K/GPa. Experimental 
values of dT c /dP in the range -0.7 to -2.0 K/GPa have 
been reported Jl5-H|. The large spread of experimen- 
tal data was, at least in part, attributed to differences in 
the samples studied [Q . Our estimate of the pressure 
dependence of T c is close to the upper limit of the exper- 
imental values. This, however, is not surprising because 
we used the large Griineisen parameter of the dominant 
E2g(T)-E2u(A) phonon branch. Non- negligible contribu- 
tions from other vibrational modes with smaller pressure 
dependences are expected to reduce the pressure shift of 
T c compared to our estimate. 

Other effects may also be of relevance, especially a pos- 
sible pressure dependence of the electron-ion matrix el- 
ement / which enters the electron-phonon coupling con- 
stant A and was assumed to be constant here. Never- 
theless, the present results (in particular the large mode 
Griineisen parameter of the phonon branch with the 
strongest electron-phonon coupling) underline the domi- 
nant role of the pressure dependence of the phonon spec- 
trum with regard to T c . 

IV. EXPERIMENTS 

A. Experimental Details 

Polycrystalline MgB2 was prepared from stoichiomet- 
ric mixtures of Mg (Johnson Matthey Inc., 99.98%) and 
natural Boron powder (~325 mesh, 99.99%, Aldrich). 
The reaction was carried out at 850°C for two days us- 
ing sealed Ta capsules under Ar atmosphere that were 
in turn encased in evacuated silica tubes. After grinding 
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under Ar, the sample was annealed at the same temper- 
ature for one day. The x-ray powder diffraction pattern 
showed a single phase sample with a = 308.63(3) pm and 
c = 352.34(3) pm. In susceptibility measurements the su- 
perconducting transition was observed at T c — 38.2(1) K. 

For the high-pressure experiments at ambient temper- 
ature small amounts of the powder samples were placed 
into gasketed diamond anvil cells (DACs). The remain- 
der of the sample cavity was filled with KC1 as a quasi- 
hydrostatic pressure medium. The sample was in direct 
contact with the diamond window through which the Ra- 
man and reflectance spectra were taken. 

High-pressure Raman spectra were excited at 633 nm 
employing a long-distance microscope objective. They 
were recorded in back-scattering geometry using a single- 
grating spectrometer with a multi-channel CCD detector 
and a holographic notch filter for suppression of the laser 
line (Dilor LabRam). For the Raman experiments the 
DAC was equipped with synthetic diamonds (Sumitomo 
type Ha) which emit only minimal luminescence. Op- 
tical reflectivity spectra of MgB2 in the energy range 
0.6-4.0 eV were measured using a micro-optical setup 
described in Ref. |47]]. Pressures were measured by the 
ruby luminescence method W8 49 1 . Low-temperature Ra- 



man spectra of MgB2 were excited with the 647-nm line 
of a Kr-ion laser. The scattered light was analyzed by a 
triple-grating spectrometer (Jobin-Yvon T64000) in com- 
bination with a multi-channel CCD detector. A DAC 
with synthetic diamonds was used as a low-luminescence 
container for the powder sample. 



B. Optical reflectance under pressure 

We first take a brief look at the optical reflectance spec- 
tra of MgB2, because they provide some information for 
understanding the evolution of the Raman spectra under 
pressure. At low pressures, a low reflectance is observed 
throughout the spectral range from 0.6 to 4 eV [Fig. |(a)]. 
At higher pressures the NIR reflectance increases and 
above 10 GPa we observe a Drude-like edge characteristic 
of a metallic sample. This edge then remains present even 
upon pressure release from 25 to 1 GPa. This observa- 
tion indicates that the grains of the MgB2 powder may be 
covered with some non-metallic surface layer which is, at 
least partially removed by the deformation of the grains 
upon compression. X-ray photoemission studies have in- 
deed shown that MgB2 samples which have been exposed 
to air are covered with surface layers of B2O3 |5Q] and 
possibly MgC0 3 or Mg(OH) 2 JSlJ . Werheit et al. con- 
cluded that samples of boron-rich compounds in general 
tend to be covered with surface contaminants which are 
not transparent near 500 nm [ j52j . Consequently Raman 
spectra of such compounds excited in the visible spec- 
tral region may be dominated by contributions from the 
surface layer rather than the bulk material. 



The energetic position of the Drude edge uncovered 
by pressure cycling is probably determined by the onset 
of strong interband absorption at energies around 2 eV. 
The fact that the edge shows a small red shift of less 
than 10 meV/GPa on releasing the pressure from 25 to 
1 GPa indicates a small change of the related interband 
absorption threshold. We have calculated the interband 
contributions s xx {oj) and e zz (u)) to the imaginary part 
of the dielectric function of MgB2 [see Fig. |^(b)] using 
the WIEN97 code |53| with parameters as described in 
Ref. |2(|. Without going into details here, we just note 
that the calculated £ xx (uj) indeed shows a threshold at 
2.1 eV, which shifts under pressure to higher energy at 
a rate of 8 meV/GPa. These interband transitions may 
result in resonance effects in Raman scattering. 



C. Raman spectra under pressure 

Figure [| shows ambient-temperature Raman spectra 
of MgB2 recorded at increasing pressures in the range 
0-14 GPa and decreasing pressures of 15-0 GPa. In be- 
tween the pressure was increased to ~20 GPa. The 0.6- 
GPa spectrum of Fig. ^|(a), consisting of a broad peak 
near 600 cm -1 with a high-energy shoulder, is represen- 
tative for the ambient-pressure Raman spectra of numer- 
ous samples we have investigated. It is also quite sim- 
ilar to an ambient-pressure Raman spectrum of MgB2 
reported by Bohnen et al. p3[ . With increasing pressure 
the shoulder gains intensity and the spectrum exhibits 
two clearly resolved peaks. The enhanced intensity of 
the higher-energy peak persisted upon pressure release 
until the cell was opened and the sample was exposed to 
air [Fig. |(b)] . 

Neither of the observed peaks can be attributed to 
the Raman-active E2 g (T) mode nor to the silent Bi g {Y) 
mode. The observed pressure-induced peak shifts [Fig. [7], 
Tab. |v) translate to mode Griineisen parameters of 
71 = 0.27-0.45 and 72 = 0.8-1.2 for peaks (1) and (2), 
respectively, using Bq = 145 GPa. Hence, the calcu- 
lated Griineisen parameter of the E2 g (T) phonon and 
the experimental value for peak (1) differ by a factor 
of ~6. Similarly, for peak (2) the mode Griineisen pa- 
rameter and the 0-GPa frequency differ by a factor 2.3 
and >200 cm" 1 , respectively, from the values calculated 
for the i?2 S (r) mode. The pressure- induced peak shifts 
appear not to be fully reversible. This may, at least in 
part, be attributed to difficulties in determining the ex- 
act peak positions because of an intense nonlinear back- 
ground that changes with pressure. It does not, however, 
compromise the above conclusion. 

Tentatively, we attribute the higher-energy mode (2) 
with u>o = 750 ± 20 cm -1 to a well-defined peak in the 
phonon density of states near 730 cm -1 according to cal- 
culations by Kong et al. j^[ and Bohnen et al. |^3| as well 
as inelastic neutron scattering data by Osborn et al. pi) . 
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Observation of this density-of-states peak in our spectra 
implies a violation of the momentum conservation for the 
Raman scattering process which may originate from dis- 
order of the sample and/or the small optical penetration 
depth of the laser light into the metallic sample. A vio- 
lation of the momentum conservation would not lead to 
spectral features near the frequency of the E2 g {T) mode, 
as there is no density-of-states peak at this frequency 

By taking into account the observations made in re- 
flectance measurements we conclude that only the higher- 
energy peak (2) in our Raman measurements origi- 
nates from MgB2 which is uncovered from surface lay- 
ers through the application of high pressure. The lower- 
energy peak (1) which is the dominant feature of the Ra- 
man spectra of MgB2 at ambient conditions is attributed 
to a surface layer of different chemical composition. It 
should be noted that we observed this mode also for sam- 
ples that had been stored under inert conditions. Prob- 
ably, the surface layers form already during the material 
synthesis. 

Goncharov et al. have also studied MgB2 at high pres- 
sures by Raman spectroscopy ]l9| . In contrast to our 
results they observed only a single broad spectral fea- 
ture (FWHM -300 cm" 1 ) centered near 620 cm 1 (at 
GPa) showing a large pressure dependence (7 = 2.9). 
They assign the feature to first-order Raman scattering, 
but also consider other possible interpretations. Their 
observations may be reconciled with ours if one assumes 
that the two components visible in our spectra (FWHM 
—180 cm -1 ) contribute to the very broad spectral fea- 
ture observed by Goncharov et al. The change in the 
relative intensities of the components (1) and (2) would 
then lead to a large pressure dependence of the "smeared- 
out" peak. The weighted average ( "center of gravity" ) of 
the two peak positions exhibits an effective Griineisen pa- 
rameter 7 = 2.0 [Fig. [jj. A direct interpolation between 
peak (1) at GPa and peak (2) at 14 GPa results in 
7 = 4.0. The resulting range for the effective 7 of 2.0-4.0 
thus includes the value 7 = 2.9 given in Ref. 19 . 

D. Low-temperature Raman spectra 

We investigated whether it is possible to observe a pair- 
breaking excitation below T c similar to what has been 
reported for conventional superconductors such as NbaSn 
and V 3 Si @. 

Upon cooling of a MgB2 sample below 50 K we ob- 
served a buildup of additional scattering intensity in the 
range 50-300 cm" 1 [Fig. ||(a)] which can be modeled by 
two Gaussian peaks located at 128 and 226 cm" 1 , re- 
spectively, in the 2-K spectrum. These peak positions 
are higher in energy than 2Ao ~ 90 cm -1 according to 
BCS theory, i.e. 2A » 3.5fc B T c with T c = 38 K. A num- 
ber of experiments, however, have indicated supercon- 



ductivity in MgB2 to be in the medium or strong-coupling 
regime with 2A /k B T c up to 5 (2A w 130 cm" 1 ) |§. 
With increasing temperature these peaks vanish only at 
a temperature between 74 and 100 K, i.e., well above T c 
[Fig. ||(b)]. They were observed only for parallel polar- 
izations of the incident and scattered light. Furthermore, 
the rising intensity below 30 cm" 1 is not due to residual 
stray light but it originates from Raman scattering by 
the sample. The physical origin of the reported spectral 
features is unclear at present. 

Chen et al. also investigated the superconducting gap 
of MgB2 by Raman spectroscopy plj . They reported 
a redistribution of spectral weight at frequencies be- 
low —200 cm" 1 and the appearance of a small peak at 
110 cm -1 which they attributed to the pair-breaking ex- 
citation. Differing from our results a polarization depen- 
dence of the Raman spectra was not observed. 

V. CONCLUSIONS 

Using the frozen phonon approach we have calculated 
the harmonic frequencies of selected optical phonons of 
MgB2 for different volumes corresponding to a pressure 
range of to 30 GPa. Large differences are found in the 
absolute values of mode Griineisen parameters. In partic- 
ular, for the whole Ei g (T)-E2 U (^4) branch the frequencies 
are highly sensitive to volume changes, as indicated by 
large values of the corresponding mode Griineisen param- 
eters. In earlier calculations ||,^3| the E2 g (T)-E2 U {A) 
branch of the phonon dispersion was identified as the 
one which exhibits the strongest electron-phonon cou- 
pling. Therefore, within a medium or strong coupling 
scenario it is considered to govern the superconducting 
properties of MgB2. Within this picture the present re- 
sults for the mode-Griineisen parameters underline the 
dominant role of the phonon spectrum with respect to 
the observed decrease of T c with pressure. 

Based on our calculations and Raman scattering exper- 
iments on MgB2 at high pressures we conclude that the 
dominant Raman peak near 600 cm -1 in the ambient- 
pressure spectra reported by various authors does not 
originate from MgB2. We rather relate it to a contami- 
nent phase present at the sample surface. A second peak 
located near 750 cm -1 , visible only as a weak shoulder 
in spectra of as-grown samples, is tentatively attributed 
to a peak in the phonon density of states. A calculation 
of the phonon density of states under pressure, extending 
previous ambient-pressure work [jslps]] , is desirable to test 
this assignment. 

Temperature-dependent Raman experiments on MgB2 
showed a build-up of additional scattering intensity in the 
region below 300 cm" 1 upon cooling below T c . A broad 
double-peak structure appeared at 130/230 cm" 1 and 
vanished only upon warming the sample to temperatures 
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well above T c . Further study of the low-temperature Ra- 
man spectra of MgB2 is clearly needed to elucidate the 
physical origin of the various observations. 
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TABLE I. Calculated structural parameters of MgB2 at zero pressure: equilibrium volume (Vb), lattice constants (ao,co), 
bulk modulus (-Bo) and its pressure derivative (B'). Variation of the calculated c/a ratio with volume is given by the quadratic 
polynomial c/a = co/fflo + ct(l — V/Vo) + — V/Vo) with the coefficients a, /? listed below. 





V 


Oo 


CO 


Bo 


B' 


c /a 


ot 







(A 3 ) 


(A) 


(A) 


(GPa) 










Present calc. 


28.663 


3.065 


3.522 


144.7 


3.6 


1.1490 


— 0.258 


-0.06 


Experiment 


28.917(f) a 


3.08230(2) a 


3.5f461(5) a 


147-155 b 


(4.0 b ) 


1.14026 a 






Difference 


-0.9% 


-0.6% 


+0.2% 












Previous calc. 


28.888 


3.075 


3.527 


140.1 


3.93 


1.1468 


— 0.211 


-0.272 



a Ref. |l|L at 37 K. 

b Refs. |^,^|Jl|] with the assumption B' = 4. 



TABLE II. Calculated pressure and volume dependence of selected phonon frequencies in MgB2. The zero pressure frequency 
loo and the linear and quadratic pressure coefficients were obtained by least square fits of lj(P) — too + a ■ P + b ■ P 2 to the 
calculated data, and the mode Griineisen parameters 70 (at equilibrium volume) are derived from a similar quadratic expression 
for to(V). P has been obtained from V through the calculated P — V relation (see text). 



Mode 


Wo 


a 


b 


7o 




(cm" 1 ) 


(cm^/CPa) 


(cm-VGPa 2 ) 






695 


3.06506 


-0.018995 


0.6 




535 


8.9744 


-0.0780 


2.5 


B 2 u(A) 


636 


1.70378 


-0.01163 


0.4 


E 2u (A) 


480 


8.9019 


-0.0809 


2.8 



TABLE III. Harmonic contributions to the phonon frequencies at the F-point calculated by different ab initio methods. All 
frequencies are given in cm- 1 . 



Bi 9 (T) 



£2 S (P) 



A 2u (P) 



E lu (T) 



Method 



Reference 



690 
679 
692 
696 
702 
702 



470 
665 
585 
536 
571 
486 a 



390 
419 
401 
394 
405 
402 



320 
328 
335 
322 
327 
328 



Frozen Phonon 
Linear Response 
Linear Response 
Linear Response 
Linear Response 
Frozen Phonon 



Kortus 
Satta 
Kong 
Bohnen 
Bohnen 
Yildirim 



I : 
II: 



at V e: 
at Vb" 



695 



535 



400 " 



333° 



Frozen Phonon 



Present work 



a Frequency based on the harmonic term of the quadratic expansion of AE to t(u) (analogous to Eq. ^|). The "shifted" frequency 
obtained by treating the anharmonicity within the Self-Consistent Harmonic Approach or based on the energy levels of an 
anharmonic oscillator amounts to 565 or 601 cm" 1 , respectively. 
b Results obtained by fitting the calculated E(u) with a quadratic polynomial. 



TABLE IV. Observed frequencies luq of Raman features of MgB 2 and their pressure coefficients a = dui/alP. The values for 
the mode Griineisen parameters 70 are based on a bulk modulus Bo = 145 GPa. 



Mode 


Wo 


a 


70 




(cm" 1 ) 


(cm-VGPa) 




Peak (1) P up 


597(1) 


1.85(15) 


0.45 


P down 


609(2) 


1.51(19) 


0.27 


up/down average 


603(6) 


1.68(17) 


0.40 


Peak (2) P up 


734(3) 


6.3(3) 


1.2 


P down 


770(8) 


4.4(8) 


0.8 


up/down average 


752(18) 


5.4(10) 


1.0 
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E2g(a) E 1u (a) E 1u (a) acoustic 




FIG. 1. Phonon displacement patterns of the zone-center (T) modes in MgB^ (the AIB2 structure, point group Dgh), in the 
order of decreasing frequencies (at zero pressure, columnwise). The labels (a), (b) denote the different patterns corresponding to 
the degenerate J5-modes. For the (optic) Aiu and E\ u patterns to be eigenmodes, the displacements u, v of the two sublattices 
have to be such that the unit cell's center of mass stays in rest, i.e. u : v = Mu g : 2Mb- The E^ g mode is Raman active, the 
A2u and Ei u infrared active, the B\ g is silent. 
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E2u( b ) E 1g( b ) E 1u(b) 




FIG. 2. Symmetry determined displacement patterns for the zone-boundary modes at the ^4-point of the hexagonal Brillouin 
zone. Only one layer of MgB2 is shown in the _B-patterns, and it is understood that equivalent atoms above and below the 
layer displayed vibrate in the opposite direction. The patterns are arranged so as to belong to the same branch of the phonon 
dispersion uij{k) as the corresponding T-modes in Fig. 1; thus e.g. the three patterns in the last column are the "end-points" 
of the acoustic branches. 
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-0.04 



-0.02 0.00 0.02 

Displacement u/a 



0.04 



FIG. 3. Calculated total energy of MgE>2 at equilibrium volume as a function of boron displacements corresponding to the 
^2g(r)-(a) and -E2 9 (r)-(b) modes of Fig. 1. Full circles and open diamonds: ab initio total energy calculations; solid and 
dashed lines: quartic polynomials, Eq. 
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FIG. 4. Calculated harmonic frequencies of the two highest modes at T and A in MgB2 as a function of relative volume; 
Vo refers to calculated static equilibrium (zero pressure). The lines represent quadratic relation fitted to the calculated points, 
and they define the mode-Griineisen parameters 70 listed in Table I. Note that 7 itself depends on pressure, its variation with 
volume is shown in the inset. The hatched areas represent the estimated frequencies of the modes E2 g (T) and E2 U {A) when 
anharmonic effects are taken into account. 
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FIG. 5. (a) Optical reflectance spectra of MgB2 at ambient temperature and pressures up to 25 GPa. Rd denotes the absolute 
reflectance at the diamond-sample surface, (b) Calculated imaginary part of the dielectric function (interband transitions) at 
GPa (solid lines) and 20 GPa (broken lines) for light polarized perpendicular (e^) and parallel (s zz ) to the c axis. 
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FIG. 6. Raman spectra of MgB2 for (a) increasing and (b) decreasing pressures (T = 300 K). In between the pressure was 
increased to ~20 GPa. After pressure release the sample was in contact with air before recording the 0-GPa spectrum in (b). 
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FIG. 7. Energies of the Raman features (1) and (2) of MgB2 as a function of pressure. Solid and open symbols refer to data 
measured at increasing and decreasing pressures, respectively. Lines represent linear relations fitted to the data. 
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FIG. 8. Raman spectra of MgB2 at ambient pressure for (a) decreasing and (b) increasing temperatures. "|| Pol." and "± 
Pol." refer to parallel and crossed polarizations of the incident and scattered light, respectively. 
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